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We propose local strategies to protect global quantum information. The protocols, which are 
quantum error correcting codes for dissipative systems, are based on environment measurements, 
direct feedback control and simple encoding of the logical qubits into physical qutrits whose decaying 
transitions are indistinguishable and equally probable. The simple addition of one extra level in the 
description of the subsystems allows for local actions to fully and deterministically protect global 
resources, such as entanglement. We present codes for both quantum jump and quantum state 
diffusion measurement strategies and test them against several sources of inefficiency. The use of 
qutrits in information protocols suggests further characterization of qutrit-qutrit disentanglement 
dynamics, which we also give together with simple local environment measurement schemes able 
to prevent distillability sudden death and even enhance entanglement in situations in which our 
feedback error correction is not possible. 
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I. INTRODUCTION 

Every physical system is embedded in a neighborhood 
of other systems, also known as its environment. Generi- 
cally, the system exchanges energy and information with 
the environment through quantum interaction processes 
which lead to system-environment correlations. If we 
cannot keep track of the environmental degrees of free- 
dom we eventually lose information about the system 
on so-called decoherence processes pQ. This is a major 
obstacle for quantum computation and quantum infor- 
mation processing, since this environmental interaction 
spoils quantum data, in the sense that it disrupts qubit 
states causing errors and degrades resources such as en- 
tanglement, necessary for some quantum algorithms. In 
the quantum scale the environment was long seen as an 
entity whose action needed to be eliminated in order to 
preserve quantum information in the system. 

As an attempt to address this issue, quantum error cor- 
recting codes were developed to suppress the degradation 
induced by the environment [2HB]. Most of these codes 
rely on redundancy and global operations, meaning they 
use many physical qubits to encode the logical qubits and 
use global operations on these several qubits to correct 
for the errors and/or detect the error syndromes. Passive 
codes, named error avoiding codes, were also developed 
and those are based on decoherence free subspaces in 
which a subspace of the system does not evolve under the 
coupling to the environment [7] . The first designed error 
correcting codes used projective measurements to detect 
the errors and unitary operations to correct them. Not 
long after the first codes, and still recently, environment 
measurements and Hamiltonian feedback were employed 
in error correction protocols [8 -15 . In the context of er- 
rors caused by spontaneous emission, hybrid codes com- 
bining passive and active codes were developed [SHU], 
in which there are clearly two different kinds of errors, 
those that are avoided by encoding the quantum data in 



a decoherence free subspace and those that need a direct 
action given that the error has occurred. Such sponta- 
neous emission correcting codes relied on the ability to 
continuously monitor the environment (global measure- 
ments on independent environments in [8] or a measure- 
ment over a global environment in [10 ), on the encoding 
of the logical qubits into entangled physical qubit states 
and also on complicated global operations on several sub- 
systems for the active error correction. It was then real- 
ized that a significant redundancy reduction was possible 
when the distance between the physical systems is much 
larger than the wavelength of the emitted excitations [11] , 
which also makes the error detection a local process. Fur- 
ther improvement on the number of required resources 
was given in [13 with the aid of a global driving Hamil- 
tonian. Therefore statistically independent environments 
probably provide a better suited scenario for less de- 
manding quantum computation. However, these propos- 
als scale badly in the sense that as the number of logical 
qubits grows one needs to collectively control the behav- 
ior of many systems through intricate global operations in 
order to protect quantum information which is encoded 
into multi-qubit entangled states. As another possible 
strategy, continuous feedback has been applied to sup- 
press decoherence due to spontaneous emission [16 , and 
feedback based strategies have also been devised to cre- 
ate and protect known entangled states [17] . However, 
all the previous efficient proposals for sustaining entan- 
glement require global operations and hence do not allow 
for distant entanglement based communication [T7j . 

Let us clarify what we mean by local and global opera- 
tions. Consider a network whose nodes are quantum sys- 
tems, and that each quantum system may be composed 
of smaller subsystems: in our case three level systems 
(or qutrits). By global operation we mean any operation 
that is not a product of single operations on the individ- 
ual subsystems, and conversely, a local operation can be 
expressed as a product of such single location operations. 
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Local operations such as single qudit operations in uni- 
son can be simple to implement, however operations such 
as measurements comprising more than 2 close-by qudit 
locations can be very hard to implement. Furthermore, 
if the subsystems are located at distant nodes global op- 
eration might be impossible to implement. Therefore, 
error correcting codes relying solely on local operations 
are much more desirable. 

In [18] the authors have proposed a local decay quantum 
error correcting code also combining active and passive 
protocols in which a logical qubit is locally encoded in a 
physical qutrit. In this code the errors are detected by 
locally measuring the independent environments and the 
system is deterministically recycled through local feed- 
back. Therefore the code is 100% local and requires nei- 
ther entangled ancillas nor global operations. Similar 
qudit encoding techniques were employed for error cor- 
recting continuous variable systems and phase dumping 
in qudit systems p~9j[20]. However, for the spontaneous 
emission codes shown here and in [18] , the minimal qudit 
dimension is d = 3, which is much less demanding and 
much simpler than the codes shown in [I9j [20] . Since 
the protocol is strictly local it can also be used to pro- 
tect known entangled states shared by spatially separated 
different parties, thus allowing for entanglement based 
quantum communication between them. 

In this paper we extend these results to the diffusion 
limit, that is, we develop a similar qutrit encoding and 
feedback error correction for the case in which the mea- 
surements performed over the reservoirs present extra 
sources of noise other than the system-environment inter- 
action (such as that introduced by a classical oscillator 
in a homodyne measurement). We also extend the anal- 
ysis previously done, testing our protocol against several 
sources of inefficiency and pointing out the most impor- 
tant ones. 

Finally, it is worth saying that this feedback error cor- 
rection is based on the indistinguishability of the decay 
channels of the qutrits in a cascade-like energy struc- 
ture. In the end of the paper, we analyze the dynamics of 
qutrit-qutrit entanglement for different energy configura- 
tions of the three-level systems and analyze the possibil- 
ity of maintaing or even enhancing entanglement through 
reservoir monitoring, also extending earlier results on the 
effects of abrupt changes and sudden death of distillabil- 
ity [21] in 3 ® 3 systems. 



II. MEASURING THE ENVIRONMENT AND 
INFORMATION FEEDBACK 

We consider a global system composed of internal 
subsystems that are weakly and smoothly coupled to 
their own local and independent reservoirs. These cou- 
plings should respect typical Markov and Born approxi- 
mations [22] [23] . We assume that the system is prepared 
in a given state (usually a pure state) and that it is not 
correlated to the environment at this stage. This can 



be expressed as |^(0)) = \S)\E), with S and E (which 
will be taken as the vacuum) designating system and en- 
vironment initial states, respectively. Correlations arise 
as a result of the evolution process that generically can 
be given by \9(t)) = ]T\ ^\Si (£)) \Ei {t)). However, in 
deriving the theoretical model one assumes that the en- 
vironment is almost unaffected by the interaction (as a 
thermal bath, in our case at zero temperature) and that 
the state presents correlations only to first order 

in the interaction coupling. Another important point is 
that the environment quickly loses memory of its past 
history as compared to the time it takes for the system 
to evolve appreciably. The interaction Hamiltonian can 
be written in the interaction picture under the rotating 
wave approximation as 



H int = J2 i [ u k h U t ) 



(1) 



in which statistically independent reservoirs are repre- 
sented by b k (t) = dubkWe-^-"*)*, with b(u) 
the environment many degrees of freedom labeled by fre- 
quencies uo and U k are system operators with character- 
istic frequencies uj k - The coupling strength is already 
included in the definition of the system operators. The 
global state vector evolution can be given by the quan- 
tum stochastic Schrodinger equation [24] 

(S) d|*(t)> = {n k dBl(t) - H\dB k {t)} |*(t)>, 

k 

(2) 

with the quantum white noise or quantum Wiener pro- 
cesses B k (t) = f*b k (s)ds. Therefore with independent 
Markovian reservoirs initially at the vacuum state we 
have the operator averages (dB k (t)dBj (£')) = 6 k iS(t — 

t f )dt, and (dBl^dB^f)) = (dsj(t)) = (dB k (t)) = 0. 
This equation is in the Stratonovich form and it is con- 
venient to convert it to the Ito form with the corre- 
sponding Ito rules dB k (t)dB\(t') = 6 k iS(t — t')dt, and 
dBl(t)dBi(t') = [dBl(t)\ 2 = [dB k (t)\ 2 = 0, yielding 



(I) d|*> 



dt + Il k dB 



Il{dB k 



(3) 

Now, if we are ignorant with respect to the environment 
and focus solely on the system properties our kno wled ge 
of the system is described by the reduced state p(t) = 
tr#{| |}, whose dynamics is given by the master 

equation 



dp = c[u\pdt = {-^{ n I n ^ p} + n fc pi4 j a, 



(4) 



in which the over- line denotes an average over the reser- 
voir. By ignoring the environment, either willingly or 
not, we also lose information on the system and this evo- 
lution leads to loss of purity and entanglement of the 
initially prepared state. 
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A. Measuring The Environment 

We may choose to perform continuous selective mea- 
surements on the environment instead of averaging over 
it. By selectively collecting the measurement results 
we acquire useful information. Take the initial state 
|^(0)) = 1*5)10), with the reservoir in the vacuum state 
(only one reservoir for simplicity) and assuming there 
is nonzero excitations in the system. If the interac- 
tion Hamiltonian H int allows for excitation exchange, 
after a very small time the global system evolves to 
|tf (dt)) = vWo)|0) + VPl|*5i)|l), with p m the prob- 
ability of one (no) excitation reaching the environment. 
This can be seen by directly iterating equation (J3|, with 
|S > = [1 - ^Udt]\S}/^ and |5i> = ILVdt\S)/y/pi, 
with pi = dt(S\U^U\S). Equation Q gives us the sys- 
tem state p(dt) = po\So)(So\ + pi|*5i)(*5i|, while mea- 
suring the environment gives us |5o(i)) with probability 
Pom. This corresponds to measuring the number opera- 
tor in the environment dA(t), with A(t) = (s)b(s)ds, 
whose eigenvalues are the number of excitations counted 
in a time interval dt, one or zero excitations being the 
only possibilities [24] [25]. Such measurement scheme 
is called a quantum jump unraveling [26 . In principle, 
there are infinitely many other ways of observing the en- 
vironment. In particular, one of them is the quadrature 
operator dO = (itf + n)dt + dB*(t)+dB(t) or homodyne 
measurement that gives rise to a diffusive unraveling [27] . 
Diffusive unravelings are suitable in optics when one com- 
bines a system output with a classical local oscillator on 
a beam splitter and measures the resulting output. They 
arise naturally in condensed matter physics, for instance 
in quantum dot monitoring with a quantum point contact 
or single electron transistor [28], and are also very useful 
for entanglement [29] and geometric phase [30] charac- 
terization in open quantum systems. 

The quantum measurement results are, ultimately, 
classical records and therefore can be expressed by clas- 
sical (instead of quantum) stochastic processes. The 
eigenvalues for the measurement operators are a Pois- 
son process dN(t) for the quantum jump (where dN(t) 
assumes values or 1) or, in the case of the diffu- 
sive unravellings, dQ(t) = (U) + U)dt + dW(t), where 
dW(t) is a Wiener process [3T] . Such processes obey 
the usual Ito rules dN k (t)dN l (t / ) = 5 kl 6(t - t')dN k (t), 
dW k (t)dWi(t') = 5 kl 5(t - t')dt [32]. Thus, the condi- 
tional evolution of the system (where tilde indicates un- 
normalized state) can be written as a function of the 
measurement records 



Y,{AMdt + B ik [n k ]d^ k ( 



(*) W), (5) 



dp(t) 



with £(t) a classical stochastic function that records the 
measurement results, and represents the information we 
acquire from observing the environment. .Af[n] and 
B^[H] are superoperators acting on the density matrix 
of the system. The operators II can now be interpreted 



as a set of detectors surrounding the system in the en- 
vironment. For a jump unraveling we have AN k [H-k]p = 

-\{U\H k ,p}, and B Nk [U k ]p = U k plll - tp , with the 

Poisson process d£ k = dN k assuming the value one 
if the kth detector is triggered and zero if not, with 
mean value dN = (U^U.)dt. For a diffusive unraveling 
^Q fc [ n fc] = £[II fc ], and B Qk [U k ]p = pU\ + U k p, with 
the kth detector recording the intensity d£ k = dQ k and 
the Wiener process assuming any real value with a zero 
mean Gaussian distribution of width dt. Every detec- 
tion apparatus suffers from imperfections, which leads to 
a non unit efficiency r] not accounted for in the above 
operators. Including inefficiencies changes the jump un- 
ravelling into ANr) k [Rk]p = £[FU]p - r] k U k pUl, with a 
consequent decrease in the probability of triggering the 
detectors reflected in the mean value dN = 77 (lit II) dt. 
Whereas, the diffusive unraveling changes as dQ k = 
r] k (Tl k + H k )dt + yfrjkdW. This approach to inefficient 
detection is based only on an information perspective, in 
the sense that it describes an ability of distinguishing the 
results that decreases with the efficiency. Realistic detec- 
tion models can be found for specific systems in [33] . 



B. Information Feedback 

Now suppose that in a practical situation we wish for 
the system to behave in a controlled way, for instance 
that its state, either known or not, to be preserved. One 
strategy for obtaining such control over the system is 
quantum feedback in which the observer acts on the sys- 
tem accordingly to the measurement outcomes [34j [35] . 
Given the outcomes d£(t — r) at the time t — r the feed- 
back apparatus acts with a Hamiltonian proportional to 
an operator F at time t. The apparatus takes the delay 
time r to process the measurement result and activate 
the feedback Hamiltonian. One straightforward way [34] 
of doing this is with a Hamiltonian linear in the mea- 



-.rig(t-T) 



With the feedback 



surement outcomes — j. dt 
acting after the measurement we have the evolution of the 
system given by (to be interpreted in the Ito calculus) 



dp(t) 



t + A{:[ll}dt + Bt:[li\dZ(t) 



■ljp(t), 

(6) 

with Tp = —i[F, p\. 

For example, in the quantum jump scenario, the feed- 
back acts unitarily and in a point like manner. For non 
unit measurement efficiency 77 and nonzero feedback time 
delay the dynamics will then be given by 



dp(t) 



dN(t-r)+Cdt+B N dN(t) }p(t)-r]Up(t)UUt. 



(7) 



Whereas, in the diffusive unraveling the dynamics will be 
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'd^(t-T) 



FIG. 1. Representation of the feedback loop. The sys- 
tem emits a signal which is recorded in the photocurrent 
I(t) = by the detector II. The measurement record 

than dictates the feedback action by the Hamiltonian F. 



corrected by inefficiencies to 



dp(t) = I -dQ(t - t) + 
V 



2r] 



dt + B Q dQ(t)jp(t). 

(8) 



since tfn oc P C RR ] P C = Pc- Therefore, m a quan- 
tum jump unraveling, one can deterministically restore 
the quantum data in the codspace with a direct feed- 
back mechanism with Hamiltonian = F^jj- and re- 
cycling operation R = e _zF , if the detection operator II 
obeys the above conditions. The feedback term in equa- 
tion ([7]) can be identified with the reversing operation by 
e T p = lZp = RpR^. 




III. LOCAL RECYCLING 
A. Quantum Jumps 

Let us start by reviewing some basic concepts of stan- 
dard quantum error correction [4]. First we choose a 
subspace, a codespace, of the system to encode the log- 
ical qubits. The system is subjected to an environment 
induced error £ that can be represented in a short time 
expansion by Sp = [t + Cdt]p = ^ Tipr\. The system 
can be recycled after errors have occurred if there exists 
a recovery operation 1Z such that 1Z£ p = p for all states 
in the codespace p = PcpPc, with Pq the projector onto 
the codespace. We are interested in errors induced by 
a Markovian environment, and for the simplest case we 
have the errors To = 1 — \T^Tldt and Ti = II \fdt. Stan- 
dard error correction should then be applied as discrete 
actions over time intervals to be made very small and it 
has to simultaneously correct for different errors T^. By 
measuring and counting the excitations in the environ- 
ment we can try to design specific recycling operations for 
each of the possible errors, moreover we can even choose 
a codespace which is decoherence free with respect to at 
least one error so we do not have to correct for all of 
them. We find a codespace which is an eigenspace of the 
r operator and we have to recycle the system only when 
a quantum jump occurs. We need to find an operation 1Z 
that reverses the quantum measurement and we choose 
to do it deterministically instead of probabilistically [36] . 

Let us then recall the conditions for deterministic re- 
versibility of a quantum jump [8j [37]. The recycling 
operation is an unitary R such that RH oc Pc (with- 
out loosing generality we neglect any effect of both the 
recycling and jump operation outside the codespace), 
therefore II oc R^Pc- Note that this implies that the 
codespace has to be an eigenspace of the Tq operator 



FIG. 2. Representation of the local recycling process. The 
logical qubits state is represented by the blue circle around 
the qutrit states. The state is preserved while the system 
does not emit an excitation. When in does the encoded state 
shifts down the qutrit states and is recycled back to its original 
subspace by the feedback action. 

Now we turn to the specific example of locally pro- 
tecting quantum information against spontaneous emis- 
sion. Suppose we have physical qutrits each with levels 
{|0), |1),|2)} with jump operators n = v / 7(|l)(2| + |0)(l|). 
In such decay process, by detecting the emitted exci- 
tations, we cannot tell which of the qutrit 's levels de- 
cayed. This is a very important condition for the pro- 
tocol to work. Another favorable condition is that the 
decay rates of transitions |2) ^ |1) and |1) \-> |0) be 
the same 7. In the long run all states are then dragged 
to the ground state. Is there a subspace in which the 
quantum jump can be reversed? The answer is affir- 
mative with the codespace being {|1), |2)}. To corrob- 
orate this one just needs to compute II^II = jPc, with 
Pc = |1)(1| + |2)(2|. The recycling operation is then 
given by R = + |0)(2| an d the feedback operator 

F = A[^(nt - n) + |0)(2| - |2)(0|] with A « 1,2092*. 
This reversing operation clearly feeds the system back 
with one excitation after it has decayed. Information is 
then preserved in the codespace while the detector is not 
triggered. As the system evolves it eventually decays and 
the qutrit state might partially or completely leave the 
codespace. Then the detector is triggered and the feed- 
back mechanism activated, and the system is recycled 
back to the original state. If the measurement is com- 
pletely efficient and the feedback time delay is null then 
the logical qubit is locally and perfectly protected from 
its coupling to the environment. The recycling process 
is illustrated in fig. Now suppose we had n qutrits, 
each with its own recycling process, then we could protect 
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any state in the codespace {|1), |2)}® n , thus protecting 
n logical qubits encoded into n individual qutrits. 



B. Quantum State Diffusion 

We can start analyzing imperfections in the recycling 
protocol by looking at a very important sort of ineffi- 
ciency that manifests itself as additional noise in the mea- 
surement records, either intentionally added or not by the 
observer, with the measurement process corresponding to 
a diffusive unraveling. We refer to the diffusive limit as an 
inefficient limit since the additional classical noise erases 
part of the information regarding the underlying interac- 
tion process between system and environment, which is 
an excitation exchange. In the jump limit we are clearly 
informed whether there has been an excitation exchange 
or not, and this is not the case in the diffusive limit. How- 
ever we may also state that we are changing the way we 
detect the error syndromes by changing the way we mea- 
sure and acquire information from the environment, and 
then we show that to reach error correction the encoding 
and the feedback Hamiltonian are also modified. There- 
fore the diffusion based recycling is different to the jump 
based and represents a different error correction code on 
its own. 

Some goals were shown to be possibly achieved with 
quantum state diffusion such as stabilization of an arbi- 
trary (but known) pure qubit state [38] and even global 
error correction [13] (relying on global encoding and 
global recovery operations). We show that it is impossi- 
ble to perfectly protect an arbitrary unknown qubit state 
encoded in a physical qubit. However, by introducing the 
qutrit encoding we can achieve error correction for a logi- 
cal qubit encoded in the qutrit. Furthermore, we globally 
protect an arbitrary unknown multi-qubit state from its 
interaction with the environment through local diffusion 
based feedback. 

In the diffusive limit the unconditioned feedback mas- 
ter equation reads (with perfect measurement and feed- 
back) 



+1. Then the appropriate feedback would be given by 
F = Y - iXS such that C[U - iF] = C[M + X(t - 5)], 
and thus the codespace is annihilated by (1 — S) and is 
preserved throughout the evolution. However, the feed- 
back Hamiltonian must be Hermitian and hence S must 
anti-commute with X, SX + XS = [13 . In the qubit 
case the only local stabilizer would be the identity opera- 
tor which does not anti-commute with any other operator 
and thus it is impossible to find the feedback Hamilto- 
nian and perfectly protect one qubit with diffusion based 
feedback. 




FIG. 3. Representation of the local recycling process. The 
logical qubit state is encoded in the blue region around the 
qutrit states. The system suffers infinitesimal errors contin- 
uously and the state is damaged in an irregular way, such 
that part of the information leaves the codespace and the re- 
maining part is also disrupted. Then the error is detected (for 
example, by homodyne measurements in optical systems) and 
the feedback mechanism is activated and the state is recycled 
back to its original subspace by the feedback action. 



In the qutrit case it is possible to find a useful stabi- 
lizer that anti-commutes with the X angular momentum 
operator. The stabilizer is then defined by 5|0) = |0), 
S\2) = |2), and S\l) = — 11), and thus the codespace is 
{|0), 1 2)} with P c = |0)(0| + |2)(2|. The recycling pro- 
cess goes as follows (see also figure [3]). In a very short 
time interval the detectors register the measurement out- 
come dQ(t) and the system state evolves accordingly as 
\S Q (t + dt)) = [1 - §ntndt + UdQ(t)] \S Q (t)), which 
damages the initial state. However, immediately after the 
detection the feedback mechanism is triggered and the 
above state is multiplied by the feedback R = e - lFd Qi t ) ^ 



dp ■■ 



[U^F + Fn, p] dt + C[U - iF\pdt, (9) thus the system actual evolution is 



and to protect an arbitrary state one needs to add a 
constant driving Hamiltonian H = — ^(U^F + FIT) to 
cancel the first term and to find an appropriate feedback 
Hamiltonian F that cancels the source of decoherence in 
the second term. The jump operator can be written as 
II = M + X + iY, with M a multiple of the identity 
and X and Y are Hermitian operators. In the particu- 
lar case of spontaneous emission we have M — and X 
and Y as the renormalized angular momentum operators 
for the physical qutrit X = & (|1)(2| + |0)(1| + h.c.) and 

Y = & (-i|l)(2| -i|0)(l| +h.c). Now suppose that the 
codespace has some stabilizer operator S, which means 
that the codespace is an eigenspace of S with eigenvalue 



\S Q (t + dt)) 



^n f n - i^(I$F + FIT) + iFU 



— )dt+(U-iF)dQ(t) 



\s Q (t)) 



-ydt 



XdQ J (1 - 5) 



\s Q (t)). 



Now, as we defined the appropriate feedback hamiltonian 
and stabilizer, it is easy to see that if the initial state 
is inside the codespace (Pc|£q(£)) = |Sq(£))) then it 
remains in the codespace |Sq(£ + dt)) = |5q(£)), since 
the term with the stabilizer S annihilates the codespace 
(t-S)\S Q ) = 2(t-P c )\S Q )=0. 
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Therefore, as in the case of quantum jump based feed- 
back, the encoding into qutrits allows for local actions to 
fully protect n-qubit states, but now with a diffusive en- 
vironment measurement and the codespace {|0), |2)}® n . 
However the diffusion based recycling can be more sensi- 
tive to inefficiencies as compared to the jump based one, 
at least in single realizations. When we consider inef- 
ficiencies, for instance in the feedback time delay, the 
trajectories become highly irregular. This behavior is a 
manifestation of the lack of information in the diffusive 
unraveling which induces degradation of entanglement. 
In fact diffusive unravelings are also know to minimize 
the average entanglement when we average over trajec- 
tories [29]. 




FIG. 4. Entanglement En of a single realization with initial 
state |12) + |21) and feedback time delay r « O.7/7. We 
also show the measurement records in both qutrits and their 
delayed version indicating the feedback action. 



C. Inefficiencies 

Now we analyze the performance of the protocol sub- 
jected to inefficiencies. Since the error detection relies 
on measuring the environment the protocol heavily de- 
pends on the measurement efficiency. The feedback time 
delay must be significantly smaller than the decay time 
otherwise the systems dissipates before the feedback ac- 
tion. The measurement time resolution dt must be much 
smaller than the decay time although it has to be grater 
than the reservoir characteristic correlation time. We fo- 
cus the analysis on the jump based code since it presents 
interesting effects on single realizations. 

It is clear that global resources such as entanglement 
can also be locally preserved by the recycling proto- 
col. Therefore, distant parties on a quantum network 
can sustain entangled states against the environmen- 
tal action and perform quantum communication proto- 
cols such as teleportation [39] and dense coding [40] . 
For such purpose using robust states would increase the 
communication efficiency. For instance, states that do 
not suffer from entanglement sudden death. W type 
states are know to be highly robust and in our recycling 
protocol they would take the forms of 1 12) + 1 21) and 
1 1 12) — |- 1 121) -|— 1 21 1) for two and three qutrits respectively. 
There is an interesting mechanism that makes such states 
robust against finite feedback time delay. Let us exam- 
ine the simplest case with an initial state 1 12) + 1 21) and 
a feedback time delay r. Suppose that an excitation is 
emitted and detected from the first qutrit and the state 
jumps to 1 02) -h 1 1 1) partially leaving the codespace. Now, 
this state evolves as 1 02) -f e _7T / 2 |ll) while the feedback 
action is delayed, and when it does act the state is put 
back into the codespace 1 12) + e _7T / 2 |21) although not 
completely restored. Now, if an excitation is detected 
from the second qutrit the state evolves and by the end 
of the delay time it is again a maximally entangled state 
e _7T / 2 [|ll) + 1 20)], and then it is finally restored by the 
feedback. Thus the undesirable evolution that takes place 
during the feedback delay on one of the qutrits can be 
reversed by the delay on another qutrit. This is the best 
possible situation given the time delay. Entanglement 



would be completely degraded if a single qutrit emitted 
two consecutive excitations within the delay time and 
this would be the worst situation. In fig. Q we show 
one particular realization presenting the effects of the 
finite feedback time delay on the entanglement of a sin- 
gle experimental run. In fig. (§ we show the average 
(over many trajectories) entanglement as a function of 
both time and feedback time delay. As shown in the fig- 




FIG. 5. Average entanglement (negativity) of the jump feed- 
back as a function of time t and feedback time delay r. 

ures the greater the time delay the less effective is the 
recycling process, although a considerable amount of en- 
tanglement can still be preserved after long times and an 
appreciable time delay. It is easy to see that states like 
1 11) + 1 22) would only be degraded by a time delay. 

The protocol is much more affected by measurement 
inefficiencies as can be seen in fig. [6j The no jump evo- 
lution starts to resemble the unconditioned master equa- 
tion evolution as the measurement efficiency decreases. It 
populates the ground state component and decreases the 
quantum coherence and entanglement in the codespace. 
As a consequence the recycling process becomes less and 
less effective as the efficiency decreases. Still an interest- 
ing effect which is what makes the protocol effective at 
small inefficiencies is that when a jump occurs (followed 
by feedback) the ground state component is eliminated 
which increases both population and coherence in the 
codespace leading to an entanglement leap (this effect 
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also played an important role in elucidating a quantum 
classical transition [41]). Thus in this process a quantum 
jump always increases entanglement, and this is shown 
in fig. [7| 




FIG. 6. Average entanglement (negativity) of the jump feed- 
back as a function of time and measurement efficiency. 




FIG. 7. Entanglement of a single realization with initial state 
1 12) + 1 21) and measurement efficiency r/ ~ 0.98. 

When the detectors have a non unit efficiency the evo- 
lution leads to loss of purity and we generally have a 
two qutrit mixed state. There is no closed expression for 
the entanglement (of formation) of such states. Here we 
choose to plot the negativity [42 as an upper bound to 
the distillable entanglement. For the sake of comparison 
we plot the negativity in all the figures. 

We have stated earlier that the indistinguishability of 
the decay channels of the qutrit is of extreme impor- 
tance for the error correction to work. Now, we re- 
lax this constraint and look at the protocol effective- 
ness as the transitions become more and more distin- 
guishable. The single detection operator II is split into 
two other detectors, U 2 = [y/a\l){2\ + y/1 - a|0)(l|] 
and III = \/7 [a/1 — a|l)(2| + ^/a|0)(l|], each one corre- 
sponding to one of the transitions. If a = 1/2 the tran- 
sitions are indistinguishable and if a = 1 (or zero) the 
transitions are perfectly distinguishable. Note that by 
changing the distinguishability we also change the over- 
all dynamics (not just the measurement process), such 
that for a ^ 1/2 we have C[U] ^ £[Hi] + £[n 2 ], however 
for a = 1/2 we have C[U] = £[ni]-bC[n 2 ]. After both de- 
tectors are triggered the same feedback action is applied. 




4 1 



FIG. 8. Average entanglement (negativity) of the jump feed- 
back as a function of time and distinguishability a. 

As shown in figure ([sj) the protocol is still highly effective 
if the transitions are slightly distinguishable. However, 
as the transitions become more distinguishable the errors 
are dominant and the system loses entanglement quickly. 




4 4 



FIG. 9. Average entanglement (negativity) of the jump feed- 
back as a function of time and the transition rate f3. 

Even though the transitions may be indistinguishable 
the transition rates may differ, which makes one of them 
more probable than the other. This is the case in har- 
monic oscillators, for instance. Generically the jump op- 
erator would be given by n = y/^[\fP\l) (2 1 + |0)(1|], yield- 
ing a probability for the transition |2) — » |1) proportional 
to /3. When the rates are unbalanced the no jump (or 
no detection) evolution disturbs the chosen codespace, 
which means that our strategy could be combined with 
a "bang-bang" error correction protocol [43j [44] to pro- 
tect the logical qubits (although we do not apply the 
"bang-bang" control). Furthermore, the jump evolution 
also disrupts the encoded information in addition to tak- 
ing it outside the codespace. As the upper transition 
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becomes more and more probable the protocol becomes 
less effective. However, for unbalanced transition rates 
corresponding to harmonic oscillators (/3 = 2) the pro- 
tocol is still able to protect a considerable amount of 
entanglement, as shown in fig. (J9|. 

Another possible imperfection in the feedback mech- 
anism is an imprecision in the implementation of the 
recycling operation, such as a fluctuating Hamiltonian 
strength as in a disorder effect. For instance the Hamilto- 
nian could be given by the ideal one with a strength which 
assumes Gaussian distributed values around the ideal 
value A each time the feedback is triggered. The disor- 
dered feedback Hamiltonian would be given by (1 + 
with S a Gaussian distributed random variable with zero 
mean. This disorder effect amounts to a non perfect recy- 
cling operation that is subjected to random imperfections 
and it is given by a one parameter family of random uni- 
tary matrices: the ideal operation times a one parame- 
ter family Re~ lF ^ . This disorder recycling process often 
does not restore the state back to the codespace since 
the disordered matrices may differ considerably (for an 
appreciable variance of the disorder distribution var{5)) 
from the ideal, thus leading to a degradation of global re- 
sources during periods in which the system ceases emis- 
sion. We show the average negativity in fig. (10). It can 



be seen that the protocol tolerates small fluctuations in 
the Hamiltonian and we believe this disorder effect to be 
not as significant as measurement inefficiency and time 
delay, although for higher disorder the effect would be 
far more damaging than a time delay since the average 
entanglement decreases approximately linearly with the 
time delay and presents a decreasing curvature with the 
disorder. 




var(S) 



0.1 4 



FIG. 10. Average negativity as a function of time and fluc- 
tuations in the Hamiltonian strength. The fluctuations are 
quantified by the variance var(S) of the Gaussian distribu- 
tion centered at the ideal Hamiltonian strength A 



IV. DISTINGUISHABLE DECAY CHANNELS 
AND QUTRIT DISENTANGLEMENT 

So far, we have shown that by including a third level 
in the usual qubit protocols we may actually prevent 
losses with, relatively robust, local strategies, which was 
not possible by working strictly with qubits. The idea 
is based on cascade decay and the indistinguishability 
of the excitations emitted by different decay channels. 
However, one could ask what happens when those chan- 
nels are completely distinguishable and therefore entan- 
glement cannot be protected by the local recycling proto- 
col. In this session we address this issue by looking at the 
dynamics of entanglement of 3 (g) 3 systems for different 
internal structures. Note that in our analysis, each sys- 
tem will always be subjected to two independent decay 
channels. As we have seen, for a decaying qubit the only 
possible structure for the decaying channel is a transition 
from the excited state to the ground state and the only 
possible behaviors of the disentanglement dynamics are 
an asymptotic or a sudden death of entanglement [45] . 
We show in this session that for qutrits, on the other 
hand, the picture is richer and may present asymptotic 
entanglement even under the action of local independent 
reservoirs. 

We study three different distinguishable decay channel 
configurations for the qutrits: 

1) The cascade or E structure with jump operators 

E l = 

2) the V structure with 

Vi = 

V 2 -- 

3) the A configuration with 

Ai = V / 7T|0)(2| 
A 2 = V^|1)<2| 



7T|0)(1| 



7il0}(l| 
7^|0)(2| 



(10) 



(11) 



(12) 



as shown in Fig. 11 

As we show bellow, different configurations will present 
different dynamics. Before we proceed, please note that, 
throughout the analysis we will always consider indepen- 
dent reservoirs of the same strength, 71 = 72, with no 
loss of generality (all the important ingredients already 
show up in this simple case) and reservoir detectors that 
are able to distinguish photons emitted by each internal 
transition hence avoiding the addition of any externally 
induced entanglement in the system. Finally, we also fo- 
cus on initial states of the form = a\00) +b\ll) + c|22) 
(where a, b and c are real numbers and a ^ b ^ c ^ for 
simplicity), which, once again, with no loss of generality, 
already present all the interesting features in qutrit-qutrit 
disentanglement . 

Previous works have introduced the possible elements 
of the dynamics of disentanglement in 3 3 systems un- 
der dissipation such as precursors of sudden death, sud- 
den death of distillability, and asymptotic decay [21] . All 
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FIG. 11. Qutrit decaying structures with one forbidden or 
low probability transition. 



these elements are present here but not necessarily for all 
the configurations. Again, as in previous cases, the pres- 
ence of those elements also depend on the initial state. 
We summarize bellow the relations between the decaying 
structure, the initial state, and the resulting general as- 
pects of the disentanglement dynamics for each internal 
structure (and show some examples in Fig. 



12): 



1. in cascade - cascade structure: (a) aymptotic decay 
for a > b > c, (b) one sudden change followed by 
asymptotic decay for b > a > c or a > c > 6, 
(c) two sudden changes and asymptotic decay for 
b > c > a, and (d) two sudden changes and sudden 
death for c > b > a; 

2. in V - V structure: (a) aymptotic decay for a > b 
and a > c, (b) one sudden change and asymptotic 
decay for a > b and a < c, and (c) two sudden 
changes and asymptotic decay for a < b and a < c; 

3. in A - A: (a) asymptotic decay with asymptotic en- 
tanglement for 2a > c and 2b > c, (b) one sudden 
change and asymptotic decay with asymptotic en- 
tanglement for 2a > c and 2b < c with c 2 < 4a6, 
(c) two sudden change and asymptotic decay for 
2a > c and 2b < c with c 2 > 4a6, and (d) two 
sudden change and sudden death for 2a < c and 
2b < c. 

We also show the negative eigenvalues of the partial 
transposition of the density matrix for some examples, 
recalling that a sudden change occurs every time one of 
the negative eigenvalues becomes non negative (fig. [13] ) . 

In all the cases we study, entanglement is measured by 
negativity which for 3 (g) 3 systems estimates the amount 
of free entanglement available in the system [46] . The 
first striking feature present in Fig. [12] is that in the 
A — A configuration, entanglement may not decay to zero 
but to a non-null value, depending on the initial state (if 
(? < Aab). This interesting property can be explained 
by the fact that there is indeed an entangled dark state 
in this configuration, composed of the superposition of 
the product of the lower states in each part, which is not 
affected by the dissipative dynamics. 

Another interesting feature is that in the V structure 
only when b or c = the system may present entangle- 
ment sudden death. In this case, we have effectively a 
two-qubit system, given that one of the upper levels of 
each qutrit is initialy not populated and, at least for zero 



cc 0.5 




FIG. 12. Examples summarizing the disentanglement dynam- 
ics for the configurations E, A, and V for different initial con- 
ditions (a, b, c). 
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Eigen-Value 1 
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FIG. 13. Examples of disentanglement dynamics showing the 
negative eigenvalues of the partial transposition of the den- 
sity matrix, such that a sudden change occurs every time a 
negative eigenvalue vanishes (or becomes positive), (a)two 
abrupt changes and sudden death in Cascade structure, and 
(b) two abrupt changes and asymptotic decay in V structure 
(a = 0.179, b = 0,2386 and c = 0,9545); (c) one abrupt 
change in Cascade structure and (d) no abrupt changes and 
asymptotic non null entanglement in A structure (a = 0.2386, 
6 = 0,9545 and c = 0,1790). 



temperature reservoirs, will remain empty for all times. 
If, on the other hand, both b and c are different from zero, 
then there will always be a residual amount of entangle- 
ment in the {|11), \22)} subspace intil all the coherence of 
the initial superposition dies, which only happens asymp- 
totically in time, hence, no entanglement sudden death. 
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Naturally, if the temperature of the reservoir is different 
from zero then the system is incoherently pumped and 
entanglement dies in finite time. 

Much in the same way as in the previous sessions, we 
can also analyze the effects of continuous monitoring of 
the reservoirs in the dynamics of the entanglement of 
the system, now as a function of the configuration of the 
levels. However, this time we only focus on no jump 
trajectories since all transition-emitted photons are dis- 
tinguishable and a single decay completely identifies the 
quantum state of the decaying qutrit and definitely kills 
entanglement. When the reservoirs are monitored and 
no decay is detected, the general evolution of the system 
will be given by 



\m) 



(13) 



with = e^/ £ |^(0)), where |^(0)) is the initial 

state and H e ff is the effective Hamiltonian (of no de- 
tection) given by H eff = -|£) n ntll n , with II n be- 
ing the jump operators. Incidently, E and V config- 
uration present the same effective Hamiltonian, conse- 
quently, the same dynamics. For the initial states under 
consideration, the non-normalized state and the negativ- 
ity are respectively given by \^(t)) = a\00) + be~ 2lt \ll) + 
ce- 2 ^\22) and 



E 



N 



(b 2 



2 )e~ 



-4ft 



{abe- 2 ^ + ace- 2 -* 1 + bce~^ t ). 

(14) 

In such no-jump trajectories there is also an effect rem- 
iniscent of the optimum singlet state conversion for 
qubits [18] . Starting with a non maximally entangled 
state there is the possibility of increasing entanglement 
before it reaches the asymptotic value. But now the sys- 
tem may attain higher values of entanglement since it is 
a qutrit-qutrit and not a qubit-qubit system (fig. ([l4|)). 

The A structure on the oder hand presents a different 
effective Hamiltonian. In this case not only it is possible 
to increase the entanglement by measurement, but there 
is also the possibility of reaching an asymptotic entangled 
state. The no-jump evolution is given by \ip(t)) = a|00) + 
6|11) + ce~ 47t |22) and the negativity is then 



En — 2 



ab + (a + b)ce 



■8-yt 



-47* 



(15) 



In the asymptotic limit the negativity is N(t — »• oo) = 
f+ft2 • Note that, this behavior does not depend on de- 
cay rates. The asymptotic entanglement in no-jump and 
unconditional dynamics occurs because the A structure 
presents a 2 02 decoherence free subspace (|QQ), |01), 1 10) 
and 1 11)), as a two-qubit system (figure ((14])). The sim- 
ple monitoring scheme is then able to prevent the abrupt 
changes, including sudden death of distillability, and in 



fact, replace them with a transient increase of entangle- 
ment followed by an asymptotic preservation of entangle- 
ment through only local operations. 




FIG. 14. The figures shows the two possible entanglement 
evolution in no jump trajectories in (a) V and Cascade con- 
figuration and (b) Lambda Configuration (a = 0:1790, b =0; 
2386 and c = 0; 9545 in trajectories whose entanglement in- 
creases; and a = 0:9545, b = 0; 2386 and c = 0; 1790 in 
trajectories whose entanglement decays). 



V. CONCLUSIONS 

We conclude by stressing that the analysis done in this 
work has identified different limitations as well as ex- 
tended the possibilities of using qutrits to encode qubits 
in dissipative single copy quantum protocols. In particu- 
lar, we have shown that in the photon detection limit the 
encoding is particularly sensitive to the efficiency of the 
measurement process, whereas quite robust against time 
delay in the feedback mechanism and very robust against 
fluctuations in the feedback interaction. We have also 
studied the limit corresponding to homodyne detection 
and have shown that the encoding in qutrits allows for 
a certain degree of protection of the initial entanglement 
which is not possible at all in this regime if the qubits are 
encoded in two-level systems. Finally, we have analyzed 
the scenario in which the decay channels are orthogonal, 
hence single clicks already kill the entanglement of the 
system. We have studied different qutrit configurations 
and we have shown that in some of them the no jump 
trajectory may still preserve entanglement even asymp- 
totically in time due to the presence of decoherence free 
subspaces in the system. 
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